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Popular Summary 

When an earthquake occurs, stress is transferred from the brittle, elastic upper crust to the more 
ductile, i.e. viscoelastic, lower crust and asthenosphere. These ductile regions flow in response to 
the changes in the ambient stress level and this flow induces postseismic (“after earthquake”) sur- 
face deformations. The surface deformations are detectable by GPS and other space geodetic 
techniques so the measurement of postseismic rebound gives insight into the rheology and geom- 
etry of subsurface geological features. 

In introductory courses on plate tectonics one is often taught that there are three types of plate 
boundaries: convergent ones where the plates approach one another, divergent ones, where they 
pull apart, and transcurrent ones where they slide past one another. In fact, plate boundaries are 
often more complex than these simple end member cases. One common type of plate boundary is 
the obliquely convergent subduction zone. Here an oceanic plate, subducts, or slides under a con- 
tinental plate at a deep sea trench. In addition the continental plate slides laterally over the sub- 
ducting plate. Instead of the relative motion being perpendicular to the trench, as it is in the case 
of the purely convergent boundary, the relative motion is at an oblique angle to the trench. 

One of the consequences of oblique convergence is that earthquakes occur not only on the 
boundary between the converging plates but also along faults lying in the crust of the overriding 
continental plate. Many of these continental seismic zones exhibit lateral motion or strike-slip 
motion that takes up the trench parallel component of the relative plate motion. The viscoelastic 
rebound to earthquakes on these crustal faults is different from viscoelastic rebound to strike-slip 
earthquakes in regions where the relative motion is purely lateral because the underlying ductile 
material is in underlain by the subducting elastic oceanic plate. This more complex geometry 
alters both the pattern of ductile flow at depth and the surface deformations. Calculations pre- 
sented in the paper show that the rates of postseismic movement can differ by 1-2 cm/yr shortly 
after a moderate to large earthquake. Thus we find that although a more complicated model is 
required to interpret postseismic rebound observations in oblique convergence zones, the mea- 
surements have the potential to reveal important information about subsurface tectonic structures. 


This paper is being submitted to the Bulletin of the Seismological Society of America 
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Abstract. According to the slip partitioning concept, the trench parallel component of relative 
plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the 
overriding continental lithosphere. The pattern of postseismic surface deformation due to vis- 
coelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault 
is modified from that predicted from the conventual elastic layer over viscoelastic halfspace 
model by the presence of the subducting slab. The predicted effects, such as a partial suppression 
of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earth- 
quake, are potentially detectable using contemporary geodetic techniques. 

Introduction 

Earthquakes transfer stress from the elastic, brittle upper crust to the more ductile lower crust 
and upper mantle. These deeper regions respond by flowing in response to the change in the devi- 
atoric stress acting on the ductile material. The ductile flow results in time-dependent deformation 
not only at depth, but also at the surface of the earth. Thus postseismic deformations depends on 
both the rheology and geometry of the relevant tectonic features. The classical, two-dimensional, 
viscoelastic rebound model, as introduced by Nur and Mavko [1974], is shown for strike slip 
faulting in Fig 1. In this model the fault, extending from an upper depth, <i/,to a lower depth, d 2 , 
and is embedded in an elastic layer of thickness, H, characterized by a shear rigidity, (X;. The elas- 
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tic layer is underlain by a viscoelastic halfspace characterized by its rigidity, [i 2 , and viscosity, r). 
Taking a coordinate system with the 1 axis aligned with fault strike, the 3 axis pointing vertical, 
and the 2 axis forming the perpendicular, a particularly simple closed-form expression for the dis- 
placement, uj, is obtained when |iy = \i 2 , dj = 0, and d 2 = D [Cohen, 1999 and references therein], 
viz: 
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Note that T(t, 1 ) = l-exp(-t/2x); the subsequent terms can be derived from the recursion relation: 
T(t,n+ 1) = T(t,n) - r' l exp(-r)/n! where r = t/2x and T = r|/(X is the Maxwell time of the viscoelastic 
material. Velocity can be derived by time differentiating the displacement equation and shear 
strain by taking the derivative with respect to x 2 . 

There have been numerous analytical and numerical extensions of the preceding formulation 
that take into account more detailed horizontal layering, lateral variations in elastic and viscoelas- 
tic parameters, and non-vertical faulting. Here we wish to focus attention on a particularly impor- 
tant issue which has received little if any attention in the past: the effect of having a subducting 
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slab lying under the crustal fault. This is a common occurrence in regions of oblique convergence 
where the trench parallel component of the relative plate motion is often accommodated by strike- 
slip faulting in the overriding continental plate [Fitch, 1972, Walcott, 1978, McCaffery, 1992], A 
good example can be found in the Wellington region of the North Island of New Zealand [e.g. 
Beanland, 1995]. Here, i.e. at the Hikurangi trough off the east coast of the North Island, the 
Pacific Plate subducts beneath the Australian Plate. The obliquity increases from north to south 
and is more than 45 degrees at the latitude of Wellington. Much of the trench parallel motion is 
taken up on the major right lateral faults of the North Island Shear Belt (NISB), notably the Wair- 
arapa and Wellington Faults, located some 140-170 km inboard of the trench. Both of these right 
lateral faults manifest long term slip rates on the order of 1 cm/yr and both are known to be seis- 
mically active with the Wairarapa Fault rupturing in a magnitude 8+ event in 1855. Thus, the pos- 
sibility of postseismic deformation, particularly viscoelastic rebound, following a strike-slip 
earthquake here, as in other zones of oblique convergence, merits consideration. In this paper I 
report the results of numerical experiments performed to elucidate the differences between the 
viscoelastic postseismic rebound patterns predicted by the simple layer over halfspace model 
(layered model) and those predicted by a model that takes the subducting plate geometry into 
account (plate model). I take the occurrence of strike-slip faulting at oblique convergence zones 
as an empirical fact. By contrast, Braun and Beaumont [1995] have previously reported the 
results of an investigation into how strain partitioning arises in elastic-brittle material driven by 
basal stresses. 

Model, Results, and Discussion 

Since we are dealing with a convergent plate boundary, it is convenient to refer to the overthrust 
block as the continental plate and the underthrust block as the oceanic plate. As in Fig. 1, the fault 
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lies at x = 200 km and the surface of the overthrust (underthrust) block is at x < 0 (x > 0). I will 
occasionally refer to the region where x < 0 (x > 0) as the continental (oceanic) plate even when 
alluding to the layered model, where no subduction occurs. The solutions for the layered model 
were obtained using the analytical expressions while the solutions for the plate model were 
obtained with finite element techniques. I verified that the finite element code reproduced the ana- 
lytical results for the layered geometry. The velocities obtained for typical calculations are shown 
in Fig. 2. In the layered case, the elastic layer was 20 km thick and the fault penetrated through 
the entire elastic layer. In the plate case, both the continental and oceanic plates had a thickness of 
20 km. The fault penetrates through the entire continental plate. The dip angle of the subducting 
plate increases with increasing distance from the plate boundary in a manner consistent with typi- 
cal subduction zone geometries. As a practical matter in generating the finite element grid, I found 
it convenient to let the subducting plate thin somewhat with depth from its unsubducted thickness 
of 20 km (see Fig. 2), but this variation in thickness does not affect the conclusions presented 
herein. The Maxwell relaxation time, T, was chosen to be 3.2 years. Fig. 2a shows the surface 
velocity one-half year after an earthquake. The deformation both landward and seaward of the 
strike-slip fault is reduced in the plate model as compared to the layered model. The difference in 
the predicted velocities derived from the models is about 1-2 cm/yr per meter of coseismic slip at 
a distance of about one fault depth from the fault. This is well within the detection capability of 
modem geodetic techniques. Thus the effects of the plate geometry should be discernible for vis- 
coelastic rebound to great, large, and even some moderate magnitude earthquakes. The physical 
basis for some of the differences in model behavior is easy to understand. Relative to the layered 
model, deformation is reduced in the plate model on the trenchward side of the fault because the 
coseismic stress change is partially supported by both the continental and underlying oceanic 
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plates. Thus the effective elastic thickness is greater than in the layered model. In fact, Fig 2b 
shows a comparison between the postseismic velocities from the plate model and a layered model 
in which the thickness of the elastic layer in the layered model has been doubled to 40 km with no 
change in the faulting depth of 20 km. Notice the similarity in the velocities on the oceanic side of 
the fault. Returning to Fig. 2a, deformation is reduced on the continental side of the fault because 
the viscoelastic material is confined to a wedge, bound on the lower side by the top of the sub- 
ducting plate. Note, however, that the fastest moving points are farther from the fault in the plate 
model than in the layered model. Fig. 3a shows the cumulative postseismic displacement after 
200 yrs by which time the displacements have diffused into the far field far more successfully in 
the layered model. At 100 km from the fault, the displacement in the layered model is more than 
80% of the coseismic slip, whereas in the plate model it is only about 50% or less. 

In Figs. 2a and 3a the thicknesses of the continental plate and oceanic plates were approximately 
equal. I now examine in more detail the effect of varying the elastic plate thickness by looking at 
end members where the elastic thickness of one of the two plates is infinite, but the fault depth 
remains 20 km. The results are shown in Fig. 3b. When the oceanic plate is thick, the postseismic 
displacements at the surface are suppressed on both the oceanic plate and the oceanward side of 
the continental plate. However, displacements on the continental side of the strike-slip fault, may 
be enhanced over that predicted from a layered halfspace model The enhancement in the near 
field is evidenced in the figure. Similarly, if the continental plate is made very thick, deformation 
on the continental side of the strike-slip fault is suppressed and deformation on the oceanward 
side of the continental plate may be enhanced or suppressed relative to the layered model depend- 
ing on the distance from the fault and the time elapsed since the earthquake. Motion at the surface 
of the oceanic plate may also be reduced compared to the layered model. 
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Many of the strike-slip faults that occur in subduction zones are not vertical structures, but, I do 
not consider this detail in the present paper for there is no simple solution comparable to equation 
(1) for non-vertical geometries. However, I do speculate that another difference between strike- 
slip faults in obliquely convergent and transcurrent zones, is that the former may extend to greater 
depths. Darby and Beanland (1992) fit the coseismic vertical defamation data (mostly beach ter- 
race heights) from the 1855 Wairarapa earthquake with a listric fault model in which slip extends 
to a depth of 40 km, much deeper than the slip in, for example, earthquakes along the transcurrent 
San Andreas Fault system. For the slip partitioning model to work in the sense of accommodating 
all the trench parallel motion on a strike-slip fault, the slip must extend to the depth at which the 
megathrust is locked. For typical subduction zones, this depth is 40 ± 5 km [Tichelaaar and Ruff, 
1993]. However, it is not necessary that the slip at all depths be coseimic. These considerations 
are relevant to the present paper in that a fault depth of 20 km that was used for the illustrative cal- 
culation is slightly deep for strike-slip faulting at transcurrent plate boundaries, but is representa- 
tive or even a bit shallow for the similar faulting at obliquely convergence boundaries. 

Summary 

Qualitatively there are many similarities between viscoelastic postseismic rebound to strike-slip 
faulting in transcurrent and obliquely convergent environments. In both cases the sense of direc- 
tion of the horizontal coseismic and postseismic motions are the same; although, this general 
behavior will be somewhat modified if the oceanic and continental plates have substantially dif- 
ferent elastic thicknesses resulting in a horizontal variation in the depth to the viscoelastic mate- 
rial. Similarly, in both cases the displacements diffuse outward with time from the fault region 
into the plate interior. However, the details of the temporal and spatial pattern of postseismic hor- 
izontal deformation are different, and well within the detection capability of contemporary geo- 



created August, 1999. Current version: September 13, 1999 1:57 pm 


7 


detic measurements for rebound from moderate to great earthquakes. Thus, the interpretation of 
postseismic rebound observations following strike-slip earthquakes in convergent zones is likely 
to require the use of numerical models that are not quite as handy as the closed-form expression 
for a layered rheology. However, properly interpreted, these geodetic measurements will have the 
potential to give additional insight into the geometry of the subducting slab and into possible dif- 
ferences between strike-slip faulting at different types of plate boundaries. 
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FIGURE CAPTIONS: 

Figure 1: Fault models. The top panel shows the layered model applicable to transcurrent plate 
boundaries and for which the postseismic displacements have been calculated using equation (1). 
The bottom panel shows the plate model applicable to oblique convergence zones and for which 
the displacements are calculated using finite element techniques. 

Figure 2: Velocity at t = 0.5 yrs versus distance from the fault for a fault slip of 1 m and a Max- 
well time of 3.2 years. The strike slip fault is located at x = -200 km. For the layered halfspace, 

H=D=20 km. (a) For the plate model, D=20 km, the dip of the subducting plate is 5.7° from the 
surface to a depth (top of the slab) of 20 km, 1 1 .3° from 20 to 40 km, 26.6° from 40 to 90 km, and 

45° from 90 to 190 km and the slab is 20 thick along a vertical path, (b) For the layered halfspace 
H - 40 km, D = 20 km. The plate model is the same as in (a) 

Figure 3: (a) Cumulative postseismic displacement after 200 yrs. (a) For the plate and layered 
models described in the text and the caption of Fig. 2a. (b) For the plate and layered models as in 
(a) and for very thick oceanic and continental plates. 
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